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ABSTRACT 

A variety of two-equation turbulence models - including several versions of the K - 1 
model as well as the K - u model - are analyzed critically for near wall turbulent flows from 
a theoretical and computational standpoint. It is shown that the K -e model has two major 
problems associated with it: the lack of natural boundary conditions for the dissipation rate 
and the appearance of higher-order correlations in the balance of terms for the dissipation 
rate at the wall. In so far as the former problem is concerned, either physically inconsistent 
boundary conditions have been used or the boundary conditions for the dissipation rate 
have been tied to higher-order derivatives of the turbulent kinetic energy which leads to 
numerical stiffness. Th e K - lj model can alleviate these problems since the asymptotic 
behavior of « is known in more detail and since its near wall balance involves only exact 
viscous terms. However, the modeled form of the u> equation that is used m the literature 
is incomplete - an exact viscous term is missing which causes the model to behave in an 
asymptotically inconsistent manner. By including this viscous term - and by introducing 
new wall damping functions with improved asymptotic behavior - a new K -r model (where 
T == 1/u, is turbulent time scale) is developed. It is demonstrated that this new model is 
computationally robust and yields improved predictions for turbulent boundary layers. 


♦ This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the author was in residence at the Institute for Computer pp ca ions in cience 
and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. INTRODUCTION 

An increasing number of practical engineering calculations of turbulent flows have been 
based on two-equation turbulence models. For many technologically important turbulent 
flows, two-equation models represent a nice compromise between zero or one equation mod- 
els and second-order closures (the former models tend to require too many ad hoc empiricisms 
whereas second-order closure models can be overly expensive for design calculations). The 
K - £ model [1] is the most popular two-equation turbulence model in use today. When 
utilized in conjunction with wall functions, the K - £ model is reasonably well-behaved and 
has been applied to the solution of a variety of engineering problems with a moderate amount 
of success. However, many important technological applications require the integration of 
turbulence models directly to a solid boundary, particularly in problems where wall trans- 
port properties are needed or where there is flow separation. The problem of developing 
low-Reynolds-number near wall corrections to the K - £ model that can be robustly and 
accurately integrated to a solid boundary remains unresolved so that models along alterna- 
tive lines continue to be proposed (see Patel, Rodi, and Scheuerer [2] for a recent review). 
Most of these near wall K - £ models involve an excessive amount of ad hoc empiricisms 
and are numerically stiff in turbulent boundary layer flows. This motivated some researchers 
to pursue alternative two-equation models based on a modeled transport equation for the 
turbulent time scale. The most notable example is the K -w model of Wilcox and co-workers 
[3, 4] where modeled transport equations for the turbulent kinetic energy K and reciprocal 
turbulent time scale u> are solved. There is considerable evidence that the K - u> model is 
more computationally robust than the K - e model for the integration of turbulent flows 
to a solid boundary. However, the K - u model yields solutions for the turbulent kinetic 
energy that are asymptotically inconsistent near a solid boundary [4]. Hence, there is the 
need to re-examine this problem from a basic theoretical and computational standpoint. 

This establishes the motivation for the present paper. 

In this paper, the near-wall asymptotics of two-equation turbulence models will be ex- 
amined from a basic theoretical standpoint. It will be shown that the K - e model has 
two major problems associated with it. The first arises from the lack of natural boundary 
conditions for the turbulent dissipation rate which has caused modelers to use a variety of 
derived boundary conditions that are either asymptotically inconsistent (e.g., the boundary 
condition of vanishing normal derivative of dissipation) or numerically stiff (e.g., the bound- 
ary condition that ties the dissipation to higher-order derivatives of the turbulent kinetic 
energy). The second problem - which can be the source of substantial inaccuracies and 
numerical stiffness - is tied to the fact that the balance of terms at the wall in the modeled 
dissipation rate transport equation depends on higher-order correlations whose models have 
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considerable uncertainties. 

It will be demonstrated that both of these problems can be largely alleviated by solving 
a modeled transport equation for the turbulent time scale r = K/e since: (a) near the wall, 
r = y 2 / 2u which provides the needed natural boundary conditions, and (b) the balance of 
terms at the wall in the modeled transport equation for r involves only exact viscous terms. 
It will be argued that these features are primarily responsible for the more computationally 
robust performance of the if -w model of Wilcox and co-workers [3, 4]. However, the K -u, 
model yields results for the turbulent kinetic energy - as well as other turbulence quantities 
- that are asymptotically inconsistent (e.g., near the wall, K ~ y 3 - 23 instead of the expected 
K ~ y2 behavior )- It will be proven that this problem arises due to the fact that an exact 
viscous cross-diffusion term is missing in the modeled w-transport equation. A new if - r 
model is obtained by including this exact viscous term and by substituting improved wall 
damping functions which are obtained by an asymptotic analysis using the results of direct 
numerical simulations of turbulent channel flow (Mansour, Kim, and Moin [5]). The new 
model will be tested for the flat plate turbulent boundary layer and comparisons will be 
made with the predictions of other models (i.e, the if - e models of Launder and Sharma 
[6] and Chien [7] as well as the if - oj model [4]) in order to assess its performance. 


2. NEAR WALL ASYMPTOTIC ANALYSIS 


For simplicity, we will restrict our analysis to incompressible turbulent flows (however, 
the crucial conclusions that will be drawn carry over to compressible flows). The mean 

velocity u and mean pressure p are solutions of the Reynolds averaged Navier-Stokes and 
continuity equations given by 


du t _ du t 
dt ^ Uj dxj 


dp 

dxi 


+ vV 2 Ui + pi 

axj 


du t 

dxi 


= 0 


( 1 ) 

( 2 ) 


where r tJ _ -u'u' is the Reynolds stress tensor, v is the kinematic viscosity, and the usual 
Einstein summation convention applies to repeated indices. We will consider the commonly 
used two-equation models based on an eddy viscosity where 


and 


S X] + ut 


( 


duA 

dxj dx i ) 


r R2 
vt = — 

€ 


( 3 ) 

( 4 ) 
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given that K s is the turbulent kinetic energy, a = pdu'Jdxjdu'Jdxj is the turbulent 

dissipation rate, and C „ is a dimensionless constant at high turbulence Reynolds numbers. 
In two-equation models, transport equations are solved for any two linearly independent 
variables constructed from K and a. In the K-t model, modeled transport equations for K 
and a are solved; in the K -w model, modeled transport equations for K and the reciproca 
turbulent time scale V s e/K are solved; and in the K — r model, modeled transport 
equations for K and the turbulent time scale r = K/a are solved. The exact transport 

equations for K and e are as follows [8]. 

(5) 

( 6 ) 


DK = T ..^I - e -V + i jV 2 K 
Dt 13 dxj 


De 

Dt 


= V e - $ e - V e + vV 2 e 


where D/Dt = d/dt + u • V. In (5) - (6), 

d (1 r-i— i . -ri' 


v =ikX ^ ww + px ) 


V, = 2v 


d ( dp 1 du 
dxi \dxj dx 


a L a r^rm 
r) + "tei 


are turbulent transport terms, and 


= 2v fciex t dx k dx t dx i dx i ozjdxjdx,. 


= 2v 


2 d 2 u[ (Pit 


(7) 


( 8 ) 


(9) 


( 10 ) 


dxjdxk dxjdxk 

are, respectively, the production and destruction of dissipation terms. 

The Taylor series expansions for the components of u { = (u ,v ,w ) are as o ows near 

wall: /auk 

u 1 = ai y + d 2 y 2 + • • • ' ' 

V 1 = biy‘ + hy 3 + ••• 

in' = Ciy + W d (^) 

where or = *(*,*,«) A = “” d “ = <*(*■*>*> giv “ that the coordinate V is normal 

to the wall (later, wall coordinates will be used where »+ = W* P™ IS the sh “ r 

vdocity). Of course, the no slip condition at the wall eliminates «o, V and c„ wh ' re as the 
continuity equation eliminates 6, (c.f, Hinse [8]). It is a straightforward matter to show a 

near a wall /, A \ 

K = 0{y 2 ), e = 0(1), r = 0{y ) v ) 
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du 

dy 


= 0(1), u' 2 = 0(y 2 ), u'v' = Q(y 3 ) 


(15) 

(16) 

(17) 

(18) 


v ' 2 = 0(y 4 ), tfS = Q(y 2 ), V = 0(y 3 ) 

V = 0(y), V c = O(y), $ e = Q(l) 

v e = 0 ( 1 ), V 2 K = 0 ( 1 ), V 2 e = 0(1) 

where 'P — r^du^/ dxj is the turbulence production. 

An asymptotic analysis of the K r model will be conducted first. In the K — e model 
the eddy viscosity near a wall is taken to be of the form ’ 

VT = (19) 

s!nlT P t‘°n C r aIySiS PreSented “ thiS Secti ° n “ dicat “ that = 0(l/») near the wall 
nee due to (15), v T must be of <% 3 ) in this region. Of course, aufflciently far from the 

wah /, assumes a value of 1. (C„ ia a conatant which ia typically taken to be 0.09) The 

p _ ( u t dK \ 

dx { l^cr*- dxi ) (20) 

where i, a conatant. From (14), (17), and (19), it ia dear that this model ia net asymp- 

tot, rally consent. However, P conaiata of two parta - the triple velocity term and the 

pressure diffusion term - aa given by (7). Direct numerical aimulationa of the Navier-Stokea 
equations indicate that 

except very close to the wall (i.e., inaide of y+ = 2; aee Manaour, Kim, and Moin [5)) and in 

th'YTT v neS lglble m COmparison t0 lhe diaaipation rate and the viacoua diffusion of 
e urbulent kinetic energy. Hence, if we approximate V by 


d . 1 . 


D = a^ ( 2 u » : > 


( 21 ) 


en the gradient tranaport model (20) ia asymptotically consistent since the right-hand-side 

o ( ) an (21) are both of 0(jr>) aa the wall is approached. Hence, it would appear that 

the asymptotic errors introduced by the use of (20) in the K — t model are probably not 
that significant. 

The turbulent transport term V. in the dissipation rate tranaport equation ia also modeled 
by a gradient transport hypothesis: 


V c = 


_d 

dx 


_ f vt de\ 

\ \ a t dx { J 


( 22 ) 
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(where cr e is a constant) in the K-e model. This model is not asymptotically consistent since 
V e = 0(1) near a wall while the right-hand-side of (21) is 0(y 2 ). However, this inconsistency 
is probably not of great consequence since both V. and $ e are of 0(1) near a wall but direct 
numerical simulations of turbulent channel flow indicate that V e « <j> t (c.f. Mansour, Kim, 


and Moin [5]). 

The production of dissipation V e and the destruction of dissipation $ e 


follows 


, £ &Ui 

V t = 


are modeled as 


(23) 


in the K - £ model where the wall damping functions /i,/a 1 awa 7 from the walL li 

is clear from (14) - (17) that these models are asymptotically consistent if /i = 0(1) and 
f 2 = Q(y 2 ) near a wall. It thus follows that the K - e model will generate solutions for 
K,e and vFv' that are asymptotically consistent if the damping functions / M = 0(1 fy) and 
f 2 = 0(y 2 ) near a wall with f\ = 1. 

While the K - £ model can be made asymptotically consistent in near wall turbulent 
flows by the introduction of only two wall damping functions - namely, = 0(1 fy) and 
f 2 = 0 (y 2) _ there are still some other major problems that need to be discussed. There 
are no natural boundary conditions on c; consequently, boundary conditions must be either 
derived or postulated. One of the commonly used derived boundary conditions is 


4 * ( 25 ) 

dy 2 

at the wall which is a rigorous consequence of the exact transport equation for K. Equation 
(25) requires information at the wall on the second-order derivative of the turbulent kinetic 
energy - a feature that can lead to considerable numerical stiffness [2]. Some of the stiffness 
can be alleviated by utilizing the alternative version of (25): 


2v 



= e 


(26) 


at the wall. However, even (26) can give rise to considerable numerical problems. 
Neumann, boundary condition 




The 

(27) 


has been used in a variety of applications of the K - £ model (c.f. Lam and Bremhorst 
[9]). Although (27) is more computationally robust, it is completely ad hoc with no solid 
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theoretical or experimental justification. In fact, recent direct numerical simulations of the 
Navier-Stokes equations for turbulent channel flow indicate that [5] 

1 de+ 

at the wall; under such circumstances the use of (27) could lead to substantial errors. 

The other major problem with the dissipation rate transport equation lies in the balance 
of terms at the wall. At a solid boundary, (6) reduces to 


dh 

"dy 2 




+ V - + Wt 


(28) 


For a fully-developed turbulent boundary layer, dcjdt = 0 and V, « as discussed earlier; 
hence, (28) simplifies to 

d 2 e , 

V = *■* ( 29 ) 

Both (28) and (29) have a major deficiency: the balance of terms at the wall involves higher- 
order correlations. This puts significant pressure on the accuracy of the near wall modeling 

of the destruction of dissipation term that can further exasperate the numerical stiffness 
problem. 

On the other hand, the turbulent time scale r = K/e, has a variety of natural boundary 
conditions. It is a simple matter to show that close to a wall 


T 



(30) 


and, hence, at the wall 

_ dr d 2 r 1 

w = V ( 31 ) 

Equations (30) - (31) have the advantage of being valid for any near wall turbulence where 
the fluctuating velocity is expandable in a Taylor series. Furthermore, the balance at the 
wall in the transport equation for r only involves exact viscous terms. This can be seen 
directly from the exact r-transport equation which takes the form 


Dt 

Dt 


T dlli T T 2 T 2 T 2 

~ 1 ~K V ~ K Vt + K*' + 


K 


, 2v dK dr 2 v dr dr 

+ Kdx t dx t ~Td^dI l +vVT 

and, hence, at the wall y = 0, the leading terms are 


(32) 


2v dK dr 2v dr dr d 2 r 

~K~dy"dy ~ T fry fry + V dtf ~ 1 = °* ( 33 ) 
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Here, each term on the left-hand-side of (33) is 0(1). The balance of terms in (33) is 
guaranteed if r = j/ 2 /2i/ near the wall. It therefore appears that the two major problems with 
the K — e model - namely, the lack of natural boundary conditions for e and the appearance 
of higher-order correlations in the balance of terms at the wall — can be overcome by the use 
of a K — r model. 

While the development of a two-equation turbulence model based on the turbulent time 
scale has been discussed in the literature (c.f. Reynolds [10] and Bardina [11]), no systematic 
study of such models has been conducted for near wall turbulence. Only the K — to model 
- which is based on a modeled transport equation for the reciprocal turbulent time scale 
u) = 1/r - has been studied in these flows to any extent (see Wilcox and Traci [3] and 
Wilcox [4]). In the K — co model, a modeled transport equation for the reciprocal time scale 
u) is solved which is of the form [4] 


Du to du x 2 , 9 

Dt ~ Cul K Tlj dx } + dxi 


(' + r) 


doj 

dxi 


(34) 


where v T = C^K/oj and C u i , C u2 and a w are constants which assume the values of 5/9, 5/6, 
and 2, respectively (again, C ^ = 0.09). However, we will now show that (34) is inconsistent 
with the exact transport equation for oj near a wall: an exact viscous term is missing and 
C W 2 must be damped. The exact transport equation for oj takes the form 


Doj V e $ e V t ojV , 2 , ojV , 2v doj dK , „ 2 

Dt K K K K ^ K ^ Kdxidxi 


(35) 


Hence, it is clear that an exact viscous cross- diffusion term - given by (2vjK )du/ dx{dKj dxi 
in (35) — is missing in the modeled cu-transport equation of Wilcox and co-workers. From 
(35), it is a simple matter to show that the leading terms in the near wall balance of oj are 


as follows: 


2 vdudK d 2 oj 2 n 

+ +oj 2 = 0 


(36) 

K dy dy dy 2 

at the plane solid boundary y = 0. Equation (36) is consistent with asymptotic solutions for 
K and oj that behave correctly near the wall, i.e., 


v • 2 ■ 2v 

K = ay , “ = ^- 


(37) 


In stark contrast to (36), the K — oj model of Wilcox based on (34) yields the balance of 
terms 

^ - C u2 oj 2 = 0 (38) 

at the wall, which is incompatible with (37). Hence, the K — oj model of Wilcox yields 
asymptotically inconsistent solutions in near wall turbulence (e.g., K ~ y 3 ' 23 ; see Wilcox 

[4])- 
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The K — u model can be made asymptotically consistent by the addition of the viscous 
cross-diffusion term 

2v dw dK 
K dxi dxi 

and by decomposition of C u 2 as follows 

C u 2 = C' u2 - 1 (39) 

where C' u3 must be damped of 0(y 2 ) near the wall. However, we feel that it is preferable 
to derive a modeled transport equation for t = l/u> since r is not singular near the wall. A 
new K — t model will be derived in the next section which is asymptotically consistent. 


3. A NEW K — t MODEL 


The exact transport equation for the turbulent time scale r takes the form 


Dt 

~Dt 


—V-1-—V-—V + — $ + —V + 
K K K t + K t + K ' + 


2v dK dr 
K dxi dxi 


2v dr dr 
t dxi dxi 


+ vV 2 t 


(40) 


which is obtained from a straightforward combination of (5.) and (6). Models for V c ,$ e , T> e 
and V are needed for closure. The production of dissipation term will be modeled as it is in 
the K — e model, i.e., 

e dui _C el dui 
V ■ = C "K Ti ’dx i = ~ Tii d^ 

where C e \ = 1.44. As mentioned earlier, this model is asymptotically consistent. The 
destruction of dissipation term, 4> e will also be modeled similarly, i.e., 



= c e 2 f 2 £ = c e 2 f 2 ^. (42) 

Here, asymptotic consistency with (17) requires that f 2 be damped of 0(y 2 ) near a wall. We 
will use a variation of the form for C e2 suggested by Hanjalic and Launder [12]: 

C e2 = 1.83[1 - jj exp(— i?e t /6) 2 ] (43) 

where Re t = K 2 / vs is the turbulence Reynolds number. Here, we set the high turbulence 
Reynolds number value of C e2 = 1.83 since it yields a somewhat more accurate value for 
the decay rate of grid turbulence than the traditional value of 1.92 (c.f. Reynolds [10]). An 
exponential form is chosen for the wall damping function f 2 as follows 


/a = [1 — exp(-y + /A 2 )] 2 


(44) 
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( 45 ) 


which is (D{jy 2 ^) near the wall. Since at the wall 


d 2 e + 

dy +2 



A 2 can be evaluated if d 2 e + /dy +2 , e + and K+ are known. By using the wall values of these 
quantities from direct numerical simulations of turbulent channel flow [5] we obtain 


A 2 = 4.9. 


The resulting model for $ e is quite similar to that proposed recently by Myong and Kasagi 

[13]. 

The turbulent diffusion term for r is defined by 


This term will be modeled by the gradient transport hypothesis 

2v T dK dr 2 vt_ dj_dr_ + _d_ f^9r\ (47) 

Vr ~ Kc t \ dx t dx x tg t2 dxi dxi dx { \a T2 dxi) 

where u T is eddy viscosity and a n and are turbulent Prandtl numbers. In deriving (47) 

it has been assumed that the turbulent transport processes parallel the molecular ones (i.e., 
each turbulent transport term is coupled with a molecular diffusion term of the same genera 
form). The turbulent Prandtl numbers for the last two terms on the right-hand-side of (47) 
are taken to be equal so that this pairing of molecular and turbulent diffusion terms is also 
true if the reciprocal time scale w is chosen as a variable instead of r (the choice of u; as a 
variable instead of r should not alter the basic physics of the model). Of course, the eddy 

viscosity is taken to be of the form 

UT = CJ^Kt ( 48 ) 

where = 0.09 and /„ is a wall damping function which is 0{ l/y) near the wall. By 
an analysis of the two distinct effects of low turbulence Reynolds number and near wall 
proximity, Myong and Kasagi [13] proposed the model 

f ll = ( 1 + 3.45//ftei)[l - exp(-y + /70)]. ( 49 ) 

This model fits the experimental data [2] reasonably well with one exception - it is asymptotes 
to one somewhat too slowly. Hence, we will consider the alternative model 

= (1 + 3.45 /yjRet) tanh(y + /70) ( 50 ) 
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since the hyperbolic tangent asymptotes to one faster by the necessary amount. 

Now, for the purposes of clarity, we will summarize the K - r model derived in this 
section: 


•ij 


2 + Vt 



VT = C^Kt 



dK dr 
dxi dx{ 


2 / i/ T \ dr dr d 

T v 77J dxidii + dx] 



(51) 

(52) 

(53) 


(54) 


0*09, £«i 1-44, while C c2> /„, and /, are given by (43), (44) and 

(50), respectively. However, to complete the model, values for the turbulent Prandtl number, 
and ”K must be provided. In this regard, we first note that if 


a n ~ °T 5 = cr K = <j e 


then the modeled r-transport equation (54) is equivalent the e transport equation 


De 


dx li 


__ — r> ' n r £ O 

Dt Gel v Tij C e2 f 2 — + — — 


K X ’d xj 


K dx. 


[(■*?)£] 


(55) 


Since the standard e-transport equation (55) is known to perform well in several high- 
Reynolds-number turbulent flows, it is reasonable to believe that <7 n ,^,a* and a e must 
assume values that are reasonably close to one another. Furthermore, for local equilibrium 
flows with zero pressure gradient and logarithmic velocity profile, we must have [2] 


Cel — C e2 — (« 2 i^e\[C~n) 


(56) 


where k — 0.4 is the Von Karman constant, 
herein, it follows that 


Hence, for the values of C el , C e2 and C M chosen 


and, hence, 


<r e = 1.36 

°n = P-m = &K = 1.36. 


(57) 
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It should be noted that the new modeled transport equation for r given by (54) is 
equivalent to the cj-transport equation 


Du) 1 X ^ du* /si r _1V) 2 

~Dt ~ ( Cel ~^K ij dxj ^ e2 ^ 2 ^ 


+2 


(-*£) 


1 OK du 


+ 


cr r i ) K dxi dxi dxi 


(v + — 

V cr T 2 


vt \ dui 
> d Xi 


(58) 


which differs from that of Wilcox and co-workers by the presence of a cross diffusion term 
and by the damping of the coefficient of u) 2 to one at the wall. 

Calculations will be presented in the next section using the common value of 1.36 for 
o n ,a n and a K which seems to be adequate for the present study. However, future research 
is needed to optimize these constants over a range of benchmark turbulent flows. 


4. COMPARISON OF THE MODELS 

Now, the performance of this new K - r model will be examined for the flat plate 
turbulent boundary layer at zero pressure gradient. Comparisons will be made initially with 
the K - w model of Wilcox [4] and the K - e model of Launder and Sharma [6] (a comparison 
with the K -e model of Chien [7] will be made later). The calculations to be presented were 
done with a two-dimensional boundary layer code based on the implicit marching scheme 
of Edwards et al. [14]. In the fully developed turbulent regime, approximately 100 grid 
points were used in the direction normal to the wall with the first grid point at y — 0.2. 
The profiles of the turbulent fields to be discussed in the figures are for a Reynolds number 
Re e « 16,000 based on the momentum thickness (this will allow for comparisons with the 
experimental data described by Patel, Rodi, and Scheuerer [2] which was compiled from a 

variety of sources including Coles [15] and Schubauer [16]). 

In Figure 1, the predictions of the K - r, K - w, and K - e model for the mean velocity 
are compared with experimental data [2], It is clear that each model yields a logarithmic 
velocity profile for 30 < y + < 300 that is well within the range of the experimental data. 
Furthermore, each model correctly yields u + = y + close to wall (i.e., for y + < 5) and 
predicts the deviations from the law of the wall for y + > 1000. In Figure 2, the Reynolds 
shear stress predicted by these three models is shown. The predictions of the various models 
are extremely close for y+ > 10. However, for y + < 10 the differences between the model 
predictions are significant. Among these models, only the K -r model yields a profile where 
vJv' ~ y 3 for y + < 10 as indicated by experiments; see Figure 3 and Patel et al. [2]. In 
Figure 4, the predictions of the K - r, K - w, and K - e models for the turbulent kinetic 
energy are compared. The K-r model yields a peak in K + of approximately 4 which is well 
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within the range of the experimental data [2] and the results of direct numerical simulations 
for turbulent channel flow [5], On the other hand, the K - u model - as well as the K - e 
model of Launder and Sharma - appear to yield peaks in the turbulent kinetic energy that 
are rather low. The turbulent kinetic energy near the wall is shown on a logarithmic plot 
in Figure 5. Only the K - r model yields K ~ t/ 2 for the entire interval 0 < y+ < 10; it 

yields the proportionality constant a+ = 0.05 - a result that is well within the range of the 
experimental data. 

In Figure 6, the profile of the turbulent dissipation rate predicted by the K — t,K — u> 
and K — e models are compared. Although the results are fairly close for y + > 20, there 
are some significant differences close to the wall. The K - t model yields a value for the 
turbulent dissipation rate at the wall of e+ = 0.1 which is quite close to the value obtained 
from experiments [2], Likewise, the peak in e + is quite close to the value obtained from 
experiments [2], In contrast to these results, the K - u, model and K - e model of Launder 
and Sharma yield values for the wall dissipation e+ that are substantially too small. In 
Figure 7, the variation of /„ with y + is shown for these three models. Only the K - r model 
is within the range of the experimental data compiled by Patel, Rodi, and Scheuerer [2], 

We did not compare directly with the results of the K -e model of Chien [7] since Chien’s 
model was calibrated by (and, hence, forced into agreement with) the experimental data for 
the flat plate turbulent boundary layer at zero pressure gradient. However, the fact that this 
model has some inconsistencies can be seen in the results for /„. In Figure 8, a comparison 
of the predictions of the K - r model and the K - e model of Chien [7] for / M is shown. It 
is clear from this figure that the model of Chien yields values for /„ that are far removed 
from the experimental data. Furthermore, the K - e model of Chien [7] has more ad hoc 
empiricisms than the K — r model presented herein. 

The skin friction predicted by the K - r model is shown as a function of the coordinate 
x along the plate in Figure 9. It is clear that the results are in excellent agreement with the 
experimental data [17], In Table 1, the fully- developed skin friction and wall dissipation rate 
are tabulated for the four models considered in this study. Only the K - r model and K - e 
model of Chien yield results that are within the range of the experimental data. However, it 

must be remembered that the Chien model was calibrated by forcing it into agreement with 
the experimental data for this flow. 

The more desirable features of r as a variable instead of e can be seen by a comparison 
of Figure 10 with Figure 6. It is clear from Figure 10 that the turbulent time scale varies 

much more smoothly with the distance from the wall; its first derivative with respect to y 
does not change sign. 
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5. CONCLUSIONS 


A basic theoretical and computational study of two-equation models for near wall tur- 
bulent flows has been conducted. The major findings of this study can be summarized as 

follows: 

(1) The K - (J model of Wilcox and co-workers [3, 4] is missing an exact viscous cross 
diffusion term. Furthermore, the destruction of dissipation term is not properly damped 
near a wall. These two inconsistencies give rise to asymptotically incorrect solutions for 
the turbulent kinetic energy (K ~ y 3 ' 23 ) near a solid boundary. 

(2) The K — £ model can be made asymptotically consistent by the satisfaction of two 
constraints: the coefficient of the destruction of dissipation term must be damped of 
0(y 2 ) near a wall, and the coefficient in the eddy viscosity must be damped of 0(1 /y) 
near a wall. Most existing corrections to the K — t model yield poor results in near wall 
turbulent flows due to the violation of these constraints. 

(3) There are numerical stiffness problems with the K — e model due to the lack of natural 
boundary conditions for the dissipation rate and the fact that the balance of terms for 
the dissipation at the wall involves unknown higher-order correlations which need to be 
modeled. These problems can, to a large extent, be overcome by the use of the turbulent 
time scale r = K/c as a variable since t has natural boundary conditions arising from 
the no-slip condition and since the wall balance for r only involves exact viscous terms. 

(4) A new K — r model was developed by making use of these ideas combined with a 
variation of wall damping functions for f 2 and /„ that were recently developed by Myong 
and Kasagi [13]. A preliminary test of this K - r model for the turbulent flat plate 
boundary layer yielded results that are quite encouraging. However, further tests and 
possible refinements are required before more definitive conclusions can be drawn. 

Future research will be directed on two fronts. The K — r model will be subjected to 
more stringent tests involving adverse pressure gradients (with possible flow separation) and 
high speed compressible flows. The turbulent Prandtl numbers a Tl and a r2 will be optimized 
over a range of such flows. In addition, the modeled r-transport equation will be used in 
conjunction with a nonlinear algebraic stress model as well as a second-order closure model. 
One of the major goals in undertaking this research was to ultimately develop second-order 
closure models that can be integrated directly to a solid boundary m complex turbulent 
flows that involve separation. These more complex issues will be the subject of our ongoing 
research effort on turbulence modeling. 
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Model 

e w 

Cf 

K-t 

0.094 

0.00245 

K - £ 

Launder- Sharma 

0.021 

0.00238 

K-t 

Chien 

0.113 

0.00244 

K -uj 

0 

0.00244 

Experiments [2,17] 

0.05-0.10 

0.00243 


Table 1. Comparison of the model predictions for the wall dissipation rate e w and skin 
friction C f ( x = 4.987m, R e = 16,465, and Re x = 1.156 x 10 7 ). 
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Figure 1. Mean velocity profile predictions for the flat plate turbulent boundary layer (Re$ ~ 
16,000): K — t model; K — u> model; • • • K — e model of Launder and Sharma 

[6]; O experimental data [2]. 
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